Scarred Resonances and Steady Probability Distribution in a Chaotic Microcavity 



in 
o 
o 

(N 



Q 

U 
d 



> 
o 
^i- 
o 
in 
o 
in 
o 



X 



Soo- Young Lee 1 ^ Jung- Wan Ryu 1:2 , Tae-Yoon Kwon 1,2 , Sunghwan Rim 1 , and Chil-Min Kim 1 

National Creative Research Initiative Center for Controlling Optical Chaos, 
Pai-Chai University, Daejeon 302-735, Korea and 
2 Department of Physics, Sogang University, Seoul 121-742, Korea 

We investigate scarred resonances of a stadium-shaped chaotic microcavity. It is shown that two 
components with different chirality of the scarring pattern are slightly rotated in opposite ways from 
the underlying unstable periodic orbit, when the incident angles of the scarring pattern are close to 
the critical angle for total internal reflection. In addition, the correspondence of emission pattern 
with the scarring pattern disappears when the incident angles are much larger than the critical 
angle. The steady probability distribution gives a consistent explanation about these interesting 
phenomena and makes it possible to expect the emission pattern in the latter case. 
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Directional lasing emission from a microcavity has 
been intensively studied in the last decade due to its 
direct application to optical communications and opto- 
electronic circuits The first realization of the direc- 
tional emission has been achieved in the microcavities 
which are slightly deformed from circle or sphere 0- In 
the slightly deformed microcavities it is known that the 
light is emitted from a boundary point of the highest cur- 
vature and the direction is tangential to the boundary. In 
ray dynamical viewpoint, this is a result of tunneling of 
the rays confined in the Kolmogorov-Arnold-Moser tori 
through the lowest dynamical barrier. 

As the microcavity is strongly deformed, the ray dy- 
namics inside becomes chaotic and lasing modes from 
the chaotic cavity can be expected to have complicated 
patterns without any directionality due to the com- 
plexity of the corresponding ray trajectories. However, 
the openness of the microcavity enhances the scarring 
phenomenon localization along unstable periodic or- 
bits, and the scarred resonance shows a strong directional 
emission from the chaotic microcavity. The scarred lasing 
modes are observed in several experimental studies 0,U 
and the emission direction from the scarred mode can be 
nontangential due to the Fresnel filtering effect 0- 

In a recent study0, it was reported that in a spiral- 
shaped chaotic microcavity there are many resonances 
showing special localized patterns, so-called quasiscarred 
resonances which have, unlike typical scarred resonances, 
no underlying unstable periodic orbit. This finding indi- 
cates that the openness of the system plays a crucial role 
in the formation of localized pattern, and implies that 
the properties of openness should be imprinted in scarred 
resonances also. It is quite interesting and important to 
uncover the mechanism of how the openness of the micro- 
cavity can change the scarred resonance pattern and can 
determine the emission direction. This is the motivation 
of our study. 

In this Letter, we illustrate characteristics of the 
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scarred resonances in a chaotic microcavity of stadium 
shape and show that the scarred resonances can be classi- 
fied into two categories, i.e., type-I and II for convenience, 
based on their difference in emission property. The 
type-I scarred resonances have a clear correspondence 
of the emission pattern with the scarring pattern, e.g., 
the bouncing point of the scarring pattern matches well 
with the emitting point. However, in the type-II scarred 
resonances there is no correspondence between emission 
and scarring patterns, and the emission pattern is deter- 
mined by the structure of steady probability distribution 
(SPD)J6| near the critical line (p c = s'mO c = 1/n, n being 
the refractive index) for total internal reflection. More- 
over, the SPD provides a consistent explanation about 
a novel feature of the type-I scarred resonances, i.e., the 
spatial splitting phenomenon of different chiral compo- 
nents of the scarring pattern. 

It is known that long time ray dynamical properties 
of a chaotic microcavity can be described by the SPD, 
P s (s,p), which is the spatial part of the exponentially 
decaying survival probability distribution in phase space 
(s,p), where s is the boundary coordinate and p = sin#, 
8 being incident angle jg. In chaotic microcavities this 
SPD is a useful tool in identifying both the structure of 
openness and the escaping mechanism of ray trajectories. 
In addition, various informations about the long time ray 
dynamics can be easily calculated from the SPD, e.g., the 
ray dynamical far field distribution is given by 

W)« f dsd P P s (s,p)T(p)6(<p- f(s,p)), (1) 

where T(p) is the transmission coefficient |7| and the far 
field angle <f> is given as f(s,p) determined by the geom- 
etry of boundary and Snell's law. This Pf ar {4>) would 
correspond to the averaged far field distribution for the 
resonances with relatively high quality factor Q. Even in 
the case of individual resonance pattern, the SPD pro- 
vides a background structure in phase space and is use- 
ful to understand the resonance pattern and its emission 
property. 

As a phase space representation of a quantum mechan- 
ical eigenfunction, the Husimi function has been widely 
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FIG. 1: (color online) The type-I scarred resonances show- 
ing a diamond (a) and a rectangle (b) scarring patterns in 
a stadium-shaped microcavity with n = Note that 

the refractive emissions come out from the bouncing points 
of the scarring patterns. The resonance positions are (a) 
kR = (91.68, -0.17) and (b) kR = (92.48, -0.035). 



used in billiard systems |2j- For dielectric microcavities, 
Hentschel et al.0 have developed the incident and the 
reflective Husimi functions as 



ff*»»c(re/Q( ) = 

Z7T 



-Fh{s,p) + {-)—ti{s,p) 



(2) 

where the weighting factor is T — \j n\p\~— p 2 and the 
components of the Husimi function are given as h(s,p) = 
fds'i,(s')Z(s';s,p), h'(a,p) = J ds'd^(s')^s';s,p), re- 
spectively. ip( s ) an( i dijj(s) are the boundary wave- 
function and its normal derivative, and £(s';s,p) is the 
minimal-uncertainty wave packet of the form 



: exp[- — - sf - ikp(s' + LI)} 



(3) 

where k and L are the wavenumber inside the cavity and 
the total length of boundary, respectively. This corre- 
sponds to a Gaussian wavefunction centered at (s,p), 
and we set the aspect ratio factor as a — {■\/2/k)~ 1 ^ 2 
throughout this Letter. This formalism for Husimi func- 
tions is exact in the high k limit, and it still gives good 
approximates around kR ~ 90 where our calculation is 
performed. 

As shown in Fig. 1, by using the boundary element 
method (BEM)[K| we obtain two examples of type-I 
scarred resonances, where scarring patterns are diamond 
(a) and rectangle (b), for a stadium-shaped microcav- 
ity with 0.2R length of linear segment, R being the ra- 
dius of two semicircles. It is clear that the emitted light 
is refracted out from the bouncing points of the scar- 
ring patterns inside the microcavity. We also obtain the 
Husimi functions H mc ( re f [ ' 1 (s,p)(E,q.$$), where the ori- 
gin of s is the right end of the stadium. Figure 2 shows 
the incident Husimi functions H mc (s,p) from which we 




FIG. 2: (color online) The H inc (s,p) of the type-I scarred 
resonances in Fig.l. (a) and (b) diamond scarring case, (c) 
and (d) rectangle scarring case. The SPD is shown by the gray 
dots and the critical line p c is denoted by the dashed line. The 
bouncing positions (s*, ±p*) expected from the diamond and 
rectangle ray periodic orbits are pointed as black dots. 



can realize two characteristics of the type-I scarred res- 
onances. First, the critical lines, p = ±p c , cross the 
localized spots, which explains why the emission pattern 
matches with the scarring pattern inside. This also im- 
plies that the type-I scarred resonances are somewhat 
leaky, i.e., relatively low Q resonances. The second is 
that peak positions (sj n± ,pj n± ) (i = 1,2,3,4) of the lo- 
calized spots in H lnc (s,p) are not located on the expected 
positions (s*,±p*) (black dots in Fig. 2) from the corre- 
sponding periodic ray orbits. They are slightly shifted in 
opposite ways depending on the chirality of propagating 
component. Moreover, the direction of the shift appears 
to be opposite in both scarred resonances, e.g., for posi- 
tive chiral component (p > 0) the shift, Ssf = s" l+ — s*, 
is negative for the diamond case and positive for the rect- 
angle case. 

In the case of a chaotic billiard, the Husimi func- 
tion of a scarred eigenfunction shows localization on 
an unstable periodic orbit and along nearby stable and 
unstable manifolds, and the resulting hyperbolic local- 
ization reflects that the system has the time-reversal 
svmmetrv|ll|. However, in the chaotic microcavity case 
the openness breaks the time-reversal symmetry and 
makes some structure in the phase space, i.e., the struc- 
ture of the SPD. The SPD structure of the present system 
is shown in the background in Fig. 2. We note that the 
expected ray positions (s*,±p*) are on the edge of the 
structure near the critical line, and the localized spots 
of H mc (s,p) are shifted toward the denser part of the 
SPD structure. This shift can be understood from the 
fact that the rays in the denser part contribute strongly 
to the interference process for formation of a resonance 
since the denser part of the SPD corresponds to the rays 
surviving longer time. 
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FIG. 3: (color online) Contour plots of H lnc (s,p) (black line) 
and H re f l (s,p)(red line) for the diamond case (a) and the 
rectangle case (b). The Goos-Hanchen shifts (c) and the an- 
gular shifts (d) with the incident angle 9i„ = arcsin \p ln± | . 
The solid diamonds and solid rectangles represent results of 
the corresponding scarred resonances. The open blue circles 
are for the hexagonally scarred resonance in Fig. 5 (a). The 
solid line is the result of a numerical calculation for a Gaussian 
incident beam at planar interface [T3|. 



A more direct explanation is possible in the H mc (s,p) 
of the rectangle case due to the simplicity of the manifold 
structure. The stable and unstable manifolds emanating 
from the rectangle periodic orbit are shown in Fig. 2 (c). 
Unlike the hyperbolic localization for a scarred eigen- 
function in billiard |TJJ, the localized spots appear to be 
deformed showing a part of the hyperbolic structure of 
the manifolds. This indicates the absence of an inflow 
along one stable manifold in the lower population region 
in the SPD due to the refractive escapes from the micro- 
cavity. As a result, the shifts to denser parts of the SPD 
not only explain the disagreement between (s*, ±f>*) and 
(s- n± ,p- n± ), but also determine the direction of the shift. 

Note that the incident angles of the scarred patterns 
are always greater than the incident angles expected from 
the corresponding periodic ray orbits, i.e., |p- n± | > p*. 
This means that the scarring pattern does not corre- 
spond to a simple connection of s* n± . In order to 
know the scar beam trajectories, we need to know po- 
sitions of localized spots in the reflective Husimi func- 
tion H re ' l (s,p). From the difference between H mc (s,p) 
and H re ' l (s,p), we can then estimate the Goos-Hanchen 
shifts |12| and angular shifts of the reflected beams, i.e., 

(Asf H± ,A P f F± ) = {\sT t± ~ af^iGC* -Pt*))- 
Since the angular shift Apf F originates from the par- 
tial refractive escape from the microcavity, it can be in- 
terpreted as the Fresnel filtering effect on the reflected 
beam and expected to be reduced as the refractive escape 
decreases. In Fig. 3 (a) and (b) are shown the differ- 
ences between H mc (s,p) and H re f l (s 1 p) for both type-I 
scarred resonances. While the Goos-Hanchen shifts are 
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FIG. 4: (color online) The positive chiral components in 
the type-I scarred resonance, (a) diamond case, (b) rect- 
angle case. The arrows correspond to the peak positions 
(s™'°"' + iP™'" 11 + ) of Husimi functions. Note that the pat- 
terns are slightly rotated from the symmetric periodic orbits. 
The negative chiral components are just the images of mirror 
reflection. 



clearly seen in both cases, the angular shift is evident 
only in the diamond case due to large refractive escape. 
The Asf H± and A8 FF are plotted as a function of in- 
cident angle in Fig. 3 (c) and (d), respectively. The 
Goos-Hanchen shifts distribute around 0.5A, A = 2ir/k, 
and the angular shifts decrease with increasing 9i n . This 
result is qualitatively consistent with the numerical cal- 
culation(solid line) for a Gaussian beam incident on a 
planar interface[13j]. 

In Fig. 4 the arrows inside the microcavity correspond 
to the peak positions {s™' ou ,p % i'' ou ) of both Husimi 
functions and the arrows for refractive escapes are ob- 
tained from the replacement of P s (s,p) by H mc (s,p) in 
Eq.JI). The background plots are the positive chiral com- 
ponents obtained by filtering the negative chiral compo- 
nents in the BEM. It is clearly shown that the inter- 
nal patterns are slightly rotated in opposite directions in 
both type-I scarred resonances, and show a good agree- 
ment with the arrows. The same plots for the negative 
chiral components are just the mirror images of the pos- 
itive ones. As a result, in the type-I scarred resonances 
the spatial splitting of two chiral components are clearly 
seen. 

While the localized spots in H mc (s,p) of type-I scarred 
resonance are located near the critical line as shown in 
Fig. 2, in type-II the localized spots are far above the 
critical line, which means that type-II scarred resonances 
have much higher Q factor due to strong confinement of 
light by total internal reflection and the refractive emis- 
sion is very weak. In Fig. 5 (a) is shown an example of 
type-II scarred resonance showing a hexagonal scarring 
pattern which is drawn on logarithmic scale for visibility 
of the weak refractive emission beams. From this figure 
we can find a striking feature of the type-II scarred reso- 
nance that light emissions do not match with the scarring 
pattern inside the microcavity. In this hexagonal exam- 
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FIG. 5: (color) (a) Type-II scarred resonance with kR = 
(91.17,-0.0137) on logarithmic scale, (b) and (c) H inc (s,p) 
on logarithmic scale and the SPD structure. The yellow dots 
represent the periodic hexagonal orbit, (s*, ±p*). (d) Far field 
distributions. The histogram is the result of the SPD and the 
red line is for the type-II scarred resonance shown in (a), and 
the blue line is for another hexagonal scarred resonance with 
higher k, kR = (170.58, -0.0117). 

pie, although the number of the bouncing points of the 
scarring pattern is six, the number of the emitting points 
is four. This mismatch between the scarring pattern in- 
side and the emitting pattern outside is very important in 
practical experiments because the type-II scarred modes 
cannot be inferred from far field measurements. 

Although the H lnc (s,p) of the hexagonal resonance 
simply shows, on the normal scale, six spots near the 



line of p = ±0.90, details near the critical line can be vis- 
ible in logarithmic scale plot shown in Fig. 5 (b) and (c). 
Here, we emphasize that the structure of H mc {s,p) be- 
low the critical line resembles closely that of SPD. This 
implies that a few of rays consisting in the hexagonal 
scarring pattern diffuses chaotically and escapes below 
the critical line and this process happens on the struc- 
ture of SPD. Note that the structure of SPD reveals the 
degree of the openness on the unstable manifold struc- 
ture. Therefore, the SPD structure near the critical line 
is the major factor to determine the emission patterns of 
type-II scarred resonances. 

The Pf ar ((f>) based on the SPD structure in Eq.Q is 
shown in Fig. 5 (d) as a histogram. This shows a good 
agreement with the far field distribution (red line) of the 
type-II resonance shown in Fig. 5 (a). We note that the 
ray dynamical analysis always corresponds to the limit of 
k — > oo. The small discrepancy of emission angles would 
be therefore reduced as k increases. The blue line in Fig. 
5 (d) is the Pf ar ((j)) of another hexagonally scarred reso- 
nance with higher k, and shows an improved agreement 
with the SPD result. 

In conclusion, we describe the properties of scarred res- 
onances in a chaotic microcavity. Their interesting be- 
haviors have been explained in terms of the SPD in ray 
dynamics. The type-I scarred resonances are character- 
ized by the correspondence between the scarring and the 
emission patterns and the spatial splitting of two chiral 
components. On the other hand, the characteristics of 
type-II scarred resonances are the mismatch of the emis- 
sion pattern with the scarring pattern inside, and the 
unique emission pattern given by the SPD. The results 
would be very crucial to understand the directionality 
and internal patterns of scarred lasing modes generated 
from a chaotic microcavity in practical experiments. 
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